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d 1 introduction 



Abstract. In this lecture I present the way stars can be modeled in two 
dimensions and especially the fluid flows that are driven by rotation. 
I discuss some of the various ways of taking into account turbulence 
and conclude this contribution by a short presentation of some of the 
first results obtained with the ESTER code on the modeling of inter- 
ferometrically observed fast rotating early-type stars. 
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Rotation remains, with magnetic fields, an ill- known quantity in the interiors of 
stars. It is however associated with fascinating objects like massive stars or the first 
stars of the Universe. These latter objects are indeed often called the factories of 
metals and because of their compactness (due to low opacity) are usually thought 
to have been fast rotators. However, understanding the mechanisms that lead to 
the enrichments of the interstellar medium with the variety of elements, requires 
the understanding of the mechanisms which mix the stars or which simply bring 
the elements from the regions of their nucleosynthesis to the surface where they 
can be expelled by winds. 

The determination of the mean flows that pervade rotating stars is there- 
fore an unavoidable step towards this understanding. This is why much work 
has been devoted to ins e rt the effects of flows into one dimensional models (e.g. 
iMaeder &: Mevne-Q 200Ct Maeder 2009t ). Even if this modeling succeeded in ex- 



plaining a variety of effects of rotation (e.g. the ratio of blue to red supergiants in 
galaxies or the abundance of lithium), ID models cannot include many specifici- 
ties of rotating fluid flows. Indeed, a fluid flow is intrinsically a multidimensional 
phenomenon. Two dimensions of space are therefore the minimum number of di- 
mensions to compute a rather general flow. For instance, the geostrophic flow of 
rotating fluids, that comes from by the domination of the Coriolis force, verifies the 
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Taylor-Proudman theorem stating that — 0. Such a condition is not compat- 
ible with the existence of a stellar core. In fact the very problem is that rotation 
imposes a cylindrical symmetry while gravity imposes spherical symmetry. The 
combination of both yields a two dimensional problem. 

In the following I therefore propose to focus our attention on the mean flows 
that pervade a rotating star. We will thus have the opportunity to go through 
the processes (baroclinicity and Reynolds stresses) that drive a secular mixing of 
the stars. Then, we will shortly present the state of the ESTER project and some 
selected first results on the modeling of early-type fast rotating stars. 



2 The problem's formulation 

In order to make a first step into the 2D-modeling of rotating stars, we shall 
consider an isolated, non-magnetic, not-mass-losing, early-type star. Moreover, we 
shall forget about any time evolution: the star is powered by nuclear reactions but 
these do not influence the chemical composition. However, there are some regions 
with turbulent flows. There we assume that the turbulence is in a statistically 
steady state. Such an ideal star does not exist but some stars like Vega may be 
close to it. 

With all the foregoing precautions, we may formulate the problem of a 2D- 
model of a rotating star in a consistent way. 

The partial differential equations that govern the steady-state of our ideal 
rotating star read: 



A(j) ^ AnGip) 

{pv ■ Vv) = -V(P) -{p)^(j) + F^ ^ ' ' 

div{pv) — 0. 

These are Poisson equation for the gravitational potential 0, the energy equation 
involving the temperature T, the density p, the entropy S, the diffusive heat flux F 
and the nuclear heat sources the momentum equation with the viscous force F^ 
and the equation of mass conservation. The brackets (...) indicate time averages. 
They are necessary to remove turbulence fluctuations. Actually, all quantities are 
time averaged. 

Equations (|2.ip should be completed by the prescriptions of the microphysics, 

namely the equation of state, the opacities, the nuclear network etc. Th ey should 

also b e supplemented with boundary conditions. All these are discussed in lEspinosa Lara fc Rieutord 
( 2013f) . Here, we shall only focus on the problems raised by the velocity and tem- 
perature fields boundary conditions. Before that we need to concentrate on the 
mean flows. 
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3 The mean flows in rotating stars 

The mean flows that pervade a rotating star and enforce some mixing have various 
origin. The most common sources are the following: 

• baroclinicity 

• Reynolds stresses 

• gravitational contraction, mass loss or mass accretion 
Let us review these phenomena in more details. 

3.1 Baroclinicity 

In fluid mechanics baroclinicity refers to the inclination of isobars and isotherms 
(or isopycnic, equipotential surfaces... ). This is indeed a natural feature of rotating 
stars that was discovered long ago bv lvon Zeipel (1924). Since much confusion has 



emerged after this work on the origin of the meridional circulation, I'd like to say 
a few words about stellar baroclinic flows in order to clarify the phenomenon and 
the way it should be approached. 

Baroclinicity emerges in fluids because pressure and temperature obey two dif- 
ferent and independent equations: there are some coupling through buoyancy and 
heat advection but these are weak. Hence, if we think to these two fields (pres- 
sure and temperature), they live independently and therefore establish there own 
system of isosurfaces. Usually, they do not coincide and baroclinicity arises. It 
arises because density depends on temperature or on both temperature and pres- 
sure. Hence, isodensity surfaces are distinct from isobars or, in other words, the 
two vectors Vp and VP are not parallel. When taking the curl of the momentum 
equation (after division by p), one gets 

^'^(^) = -2 

where the right hand side is usually called the baroclinic torque and M is a differ- 
ential operator. The consequence of this torque is that a flow arises. In stars this 
flow is basically a differential rotation, namely a purely azimuthal velocity field. 
However, as viscosity is always present, a differential rotation, like any shear flow, 
transport momentum in the direction of the velocity gradients. In a steady state 
the momentum flux induced by viscosity must be compensated by some merid- 
ional flow. Such a balance is illustrated by the ip component of the momentum 
equation, which may be written 

or, in a more condensed form, 



(3.1) 
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where we used the cylindrical coordinates (s, (p, z) and where A' represents some 
Laplacian-like operator. 

As we may note on p.ip . advection of angular momentum v ■ V(s'^f2) just 
compensates the viscous force due to angular velocity gradients I'V • (s^Vri). In 
axisymmetric situations like that of our rotating stars, advection is that of merid- 
ional circulation. We therefore see that in a steady state meridional circulation is 
controlled only by viscosity. 

Confusion arose in the past because it was assumed that rotation generates a 
thermal disequilibrium. Indeed, in a barotropic fluid 

div(xVr) ^ 

known as von Zeipel's paradox. The paradox was solved by saying that the merid- 
ional circulation velocity v was such that 

pTv ■ VS* = div(xVr) and div{pv) = 

namely, heat advection by meridional circulation compensates the thermal imbal- 
ance and verifies mass conservation. Such a reasoning is incorrect because the 
driving of a flow is due to forces or torques, not to mass and energy co nserva- 
tion w hich cannot ensure angular momentum conservation. As shown by iBusse 



()1981l) . a thermal imbalance leads to a time-evolution of the temperature field 
itself, which is more rapid than a mechan i cal rea rrangement of the fiuid. The 
foregoing discussion is detailed in iRieutordI 



3.2 Reynolds stresses 
3.2.1 The scale of turbulence 

Baroclinic flows therefore slowly advect heat, chemicals and angular momentum. 
However, the differential rotation of baroclinic origin is a mere shear flow that may 
develop instabilities if the Reynolds number is large enough or if the Richardson 
number is low enough. These instabilities usually lead to turbulence. Unlike 
thermal convection where the linear instability drives motion at scales on the 
fraction of the radius, shear instabilities in a stably stratifled fiuid inject energy 
at small-scale. To see that, we recall that the Richards on criterion taking into 



account the high thermal diffusivity of the fiuid imposes (jZahnlll992l ) 



< 1/4 

where N is the Brunt- Vaisala frequency, Sh the local shear, v and £ the velocities 
and length s cale of the eddies and k the heat diffusivity. This criterion originally 
proposed by iTownsendl (|T958), says that shear instability grows up to a scale (of 



the eddies) where heat diffusion is unable to smooth out the stable buoyancy of 
the background. Since the largest eddies have a turn-over time scale l/Sh, we may 
observe that all the eddies forced by the background shear are such that 
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T 



with 



or 



( ^ f uSh 



1/2 



where we i ntroduced R the radius of the sta r. Let us put numbers in this expres- 
sion. From lEspinosa Lara fc RieutordI (|2013f ) , we evaluate the associated turbulent 
kinematic viscosity 



ft 



SI- 
127V2 



10^ cmVs 



The global shear of differential rotation is of a fraction of the rotation rate; let us 
write Sh — f^- Hence we find 



R^ \AfR^n 



1/2 



E 



where E is the Ekman nu mber based on the shear induce d turbulence. Since 
E < 10^^ and / ^ 0.05 (see lEspinosa Lara fc Rieutordll2013l ). we find that 

i/R<3x IQ-^ 

Hence, turbulence in vertically stratified region remains on relatively small-scales 
in the vertical direction. In the horizontal directions of course the scales can be 
large (but limited by the stability of the eddies). 

As far as convection zones are concerned, nothing prevents the instability from 
driving all the scales where the Schwarzschild criterion predicts convection. Hence 
these are fully mixed. 



3.2.2 Reynolds stresses 

The convection zones are fully mixed regions, but because of the strong turbulence 
they harbour, Reynolds stresses drive mean flows like the differential rotation of 
the Sun. These stresses and the associated flows also force some flows in the 
neighbouring radiative zones. The picture that emerges from these remarks is 
that turbulence is everywhere in stars, with variable intensity, properties, and a 
modeling of its effects is required. Here we first concentrate on its momentum 
transport that is on Reynolds stresses. If we use Reynolds decomposition, namely 

V ~ (v) + if 

where v' are the velocity fluctuations with respect to the average, the Reynolds 
stress tensor is 
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R = {pv' ® v') 

Usually, correlations with density are not important (a fortiori in an incompressible 
fluid!) and the components of the Reynolds stress tensor are often written 

In a steady state, and neglecting correlations with density, mean flows therefore 
verify: 

p{v,}d,{v,} + dM<^',)) - -^^{P) + P9. + 

where F is the viscous force. 

The main question with the previous equation is the expression of the Reynolds 
stress tensor as a function of the mean-field, the so-called closure problem of mean 
field equations. 

A basic way of solving this problem is to assume that turbulent stresses are 
like fluid stresses and that they can be represented by a turbulent viscosity. We 
thus may write 



p{v,)d,{v.,)^-d.,{P)+pg^ + F2 



with 



Av + (w -v^ + 2 (yinpt ■ v 
-KV ln//t X (V X w) - ^ (^V • i/) V Inpt 



(3.2) 



where fit is the turbulent dynamical viscosity of the gas. This expression enables 
local variations of the viscosity although it does not say anything on its determi- 
nation. But there is a more fundamental difficulty; there is no reason why the 
functional form of the Reynolds stress should be 



{P'"Wj) = -Pt{di{vj) + dj{vi) - -dk{vk}S^.j) 



(3.3) 



as expected for a newtonian fluid. The analogy of turbulence with such a fluid 
is very limited because this expression is derived from the assumption that fluid 
flows are slight perturbations of the thermodynamic equilibrium. It is unclear 
what is the statistical equilibrium of turbulence and how actual flows may slightly 
deviates from it. Moreover, turbulence particles are vortices that support long 
range interactions, not collisions... 

Now, if we still admit p.3p . the expression of /if is still a problem. In the 
nineteentwenties, Prandtl introduced the idea of the mixing length suggesting that 
in turbulent shear flows 
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dvx 
dz 



where A is the mixing length and z the direction of the shear. This assumption 
leads to interesting results on wall-turbulence, but misses the point on the axis of 
a turbulent jet (where it predicts zero- viscosity!). 

Another recipe is the Smagorinski's prescription saying that fit — P^'^\{cij)\ 



where 



is the mean shear. Such a prescription is more general and still often 



used in Large-Eddy Simulations (LES) as a basic subgrid model. However, as 
the mixing-length model, this is a local prescription that does not take care of 
non-local effects. 

To avoid this pitfall, more sophisticated models have been developed for engi- 
neerin g problems. One of the most popular is the K—e model of lLaunder &: Spaldinsj 
(jl972l) : this models adds two new equations for the evolution of K the turbulent 
kinetic energy and e the energy dissipation. The advection diffusion equation that 
controls the evolution of these quantities are derived from the evolution of second- 
order correlations with a simple model of the third-order correlation (j Rieutord 
19971) . The point is that when K and e are known, the turbulent pressure and 
the turbulent viscosity are also known. These models may be useful to determine 
statistically steady inhomogeneous flows that we find in stars. However, because 
of the huge increase of the computing power, these models have been forsaken to 
the profit of LES, which are more flexible for industrial flows. In astrophysics, 
this raises the problem of the comparison with observations, which cannot be as 
detailed as laboratory experiments. 

Another line of research th at should be nientioned, is the one fo llowed by 
Riidiger and collaborators (see Riidiger 1989t Kitchatinov et al. 1994 ). This is 
the mean- field approach of turbulence where one seeks for an expression of {v[v'j) 
and other second order correlations, as functions of large-scale quantities like local 
rotation rate. In this approach the Reynolds tensor is not reduced to that of the 
functional form of a newtonian fluid. New effects like the A-effect or anisotropic 
turbulent v iscosities, appear. T hese approaches of course raise new kinds of prob- 



lems (see ISnellman et al.ll2012l) 



3.3 Gravitational contraction, mass accretion and mass-loss 

The last source of large-scale motions in a rotating star is due to expansion or con- 
traction of the star at various phase of its evolution. In non-rotating stars these 
phenomena just lead to radial flows, but in a rotating star, angular momentum 
conservation leads to a differential rotation and an associated meridional circula- 
tion. Not much is known on these flows that are currently under investigations 
(Hypolite & Rieutord in preparation). Again, the associated shear will lead to 
some small-scale turbulence in radiative regions. 
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3.4 Conclusion 

To conclude the foregoing discussion, let us underline that we considered here only 
non-magnetic processes. Magnetic fields will of course complicate the picture. An 
open question is how all these sources compete together and of course how we 
should model the associated turbulence. 



4 Heat transport 

4.1 convection zones 

The foregoing discussion focused on turbulent viscosity and more generally on 
momentum transport. However, heat transport was really the first obstacle to 
circumvent when first stellar models were constructed. The well-known mixing- 
length theory, based on Prandtl's idea, offers a plausible answer to this problem. 
In this approach, all the difficulties of turbulent heat transport are condensed in 
a dimensionless parameter of order unity q;mlt that has been calibrated on the 
Sun. Unfortunately, such a parameter is not universal and varies from star to star, 
impe ding any precise predicti on. A clear example is given by the two stars of a 



Cen ( Eggenberger et al.ll2004 ) 



In two dimensions, the difficulty is squared. First because the ID-MLT has no 
straightforward generalisation to 2D. Turbulent transport is now also in latitude 
via turbulent diffusion and mean fiows. 

A first idea is use the "down-gradient" prescription for the heat flux (this 
is a prescription that is largely used in the K — e model). Here, since thermal 
convection disappears when the entropy gradient vanishes or becomes positive (in 
the direction opposite to effective gravity) , a natural way to model the convective 
heat flux is to assume that it is proportional to the entropy gradient, namely 

Fconw = (pCpTv) = -XturhTVS/Tl 

where TZ is the ideal gas constant. 

When reduced to ID this prescription can be related to the classical MLT. 
It may therefore be viewed as a generalization of this theory. However, there 
are other jpos sibilities like t he introduction of an eddy conductivity Xt such that 
F = -xt^T (iRudigeilllQSgh . 



4.2 The surface tieat transfer 

One of the new problem raised by 2D-models concerns the implementation of 
the boundary conditions that should be met by the temperature field. A simple 
prescription that replaces a detailed atmospheric model is to impose that the star 
radiates locally as a black body, namely that 



— Xrn • Vr — aT^ 



(4.1) 
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where Xr is the radiative conductivity, cr Stefan-Boltzmann constant and n the 
outer normal of the surface. In non-rotating models this conditi on is usuall y 
applied where the optical depth is Ts = 2/3, but = 1 is also used ( Morellll997h . 



In 2D-models, the situation is more complicated: one has to define the surface 
where boundary conditions (for velocity, gravitational potential and temperature) 
are imposed. Since the opacity is a rapidly varying quantity at the surface of a 
star an iso-optical-depth surface is not convenient numerically. We therefore chose 
to impose boundary conditions on an isobar and select the isobar that "emerges" 
at the pole, i.e. whose pressure is 

P. = r.^, (4.2) 

^polc 

On this isobar T ~ Tcs at the pole only. As one moves towards the equator this 
isobar sinks into the optically thick matter because at Ts =2/3 pressure is stronger 
at the pole than at the equator (recall that gravity is larger at the pole) . The trick 
here is to determine the temperature on this isobar as a function of the latitude. 
For this we assume that the matter lying above this iso bar behaves like a polytrope 



of ind ex n. With this hypothesis, it can be shown (cf lEspinosa Lara fc Rieutord 



20131) that the temperature on this isobar follows the law 

I ffpoic \ ' I -Xrn-VT 

^9cS{e) Kpoic, 



n{9) = ^f^^ ] I : ) (4.3) 



which fully determines the temperature field. In more sophisticated models, Ti,{0) 
may be determined by the modeling of the stellar atmosphere. 



5 Computing the baroclinic flows in 2D-models 

We focus on early-type stars so as to avoid the modeling of outer convective layers. 
Thermal convection is only modeled in the core of the star. One-dimensional 
models tell us that i t is very efficient and therefore the assumption of an isentropic 
core is quite good ( Maedeil [ioog) . ESTER models therefore assume isentropic 
cores and no further modeling of core turbulence is included yet. 



5.1 Flows in the radiative envelope 

Above the core, the radiative envelope is not in hydrostatic equilibrium as re- 
called above. Baroclinic flows pervade it and are responsible of the mixing. They 
represent the first difficulty for the computation of self-consistent models. When 
determined, we have a prediction of the differential rotation and the meridional 
circulation. 

As we noted in sect. 13.11 the determination of these flows rests on the fluid's 
viscosity, which is crucial for stationary solutions. The importance of viscosity is 
shown by the value of the Ekman number, namely 
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As shown in lEspinosa Lara fc RieutordI (|2013l ). its value is always less than 10^^ 
even if shear turbulence is accounted for. 

Viscous effects are therefore expected to be small. But they are crucial as they 
lift the degeneracy of the differential rotation as we shall show now. 

Let us take the curl of the momentum equation and let us project it in the 
azimuthal direction on e^. If viscous or Reynolds stresses are neglected then one 
finds: 



VP X Vp 



(5.1) 



dz 

The solution of (|5.ip is invariant with respect to the addition of an arbitrary 
function of s; if f2(s, z) is a solution 

VL'^{s,z) = n^{s,z)+F{s) 

is also a solution. As shown in iRieutordI (|2006al ). viscosity lifts this kind of de- 
generacy. This may be understood in the following way. The differential rotation 
depends on F{s). The balance of angular momentum advection and diffusion is 
given by p.ip which is completed by mass conservation. Once $7 is known, these 
two equations give the meridional circulation, but without matching the boundary 
conditions. In general, the derived meridian flow does not verify 



u • r? = 

on the stellar surface. This condition is enforced by an Ekman boundary layer 
where the mass flux into the surface is identified with the mass pumping of the 
boundary layer. Indeed, the Ekman boundary layer completes an inviscid solution 
like sri(s, z) + G{s) in such a way that stress-free conditions are met, but the 
boundary layer correction ue^Utp usually do not verify mass conservation. This 
latter constraint is corrected by the velocity component orthogonal to the layer, 
which is called the pumping. The identification of the pumping of the layer with 
the mass flux generated by the circulation at the stellar surface gives the differential 
equation for the unknown geostrophic flow that appears in invisc id solutions {F{s) 
or G{s)). This differential equation is derived in lRieutordI (|2006al) . in the simplified 
case of an incompressible "star" not far from rigid rotation. 

In modeling rotating stars in 2D we are not interested in the boundary layer 
because they are ve ry thin and would require a high spatial resolution near the 
surface. This is whv lEspinosa Lara fc Rieutord ( 2013 ) derived a special boundary 
condition that completes (15. ip and insure stress-free conditions without an explicit 
computation of the boundary layer flow. This boundary condition reads 



E^s^i ■Vn + il)f ■ V(s^O) = at the surface . 



(5.2) 
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Fig. 1. Meridional view of the differential rotation of 30 Mq ZAMS star rotating at 98% 
of its critical angular velocity with X=0.7 and Z=0.02. 

where Eg is the surface Ekman number, ^ and f unit vectors normal and tangential 
to the surface, while is the stream function of the meridional circulation. It turns 
out that this circulation scales as the interior Ekman number, therefore this small 
parameter drops out of the momentum equation and only order one quantities are 
computed. 




5.2 The Core-Envelope Interface 

The foregoing simplification unfor tunately breaks down at the core -envelope in- 
terface (CEI). Indeed, as shown bv lEspinosa Lara fc RieutordI ( 2013 ) the chemical 



evolution of the convective core, that leads to a density discontinuity at the CEI, 
also leads to an angular velocity discontinuity on this interface. In addition we 
may also consider that turbulent viscosity steeply increases when one goes from the 
envelope to the core. These discontinuities are triggering the so-called Stewartson 
layer that develops along the tangential cylinder circumventing the core. This layer 
may play a crucial role in the chemical enrichment of the envelope by products of 
the nuclear reactions. Its computation is however demanding as it scales like E^/^ 
but fortunately less than the Ekman layer whose thickness scales like E^/^. More 
studies are necessary to better determine the properties of this interfacial region, 
which is also suspected of harbouring some convective overshooting. 



6 Some results of the ESTER project 



A detailed account of the present achievements of the E ST ER project may be found 
in the two papers bv lEspinosa Lara fc RieutordI ( 2013 ) and Rieutord fc Espinosa Lara 
(|2012l ). Here, we wish to give a brief summary of the performance and first results 
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Fig. 2. Variations of the flux with colatitude for a ZAMS star of 30 M0 rotating at 98% 
of its critical angular velocity. Solid line: the ratio of the flux to the polar flux as a 
function of the colatitude. Dashed line: the prediction of the von Zeipel law. 

of the ESTER codfl 
6.1 Presentation 

The ESTER code is solving the equation of stellar structure in two dimensions 
assuming an isentropic convective core and a radiative envelope, thus re- 
stricted to the modeling of an early-type star. We use OPAL tables for the equation 
of state and opacities. Convection in layers is not computed yet : the temperature 
gradient is assumed to be close to the radiative one, which is fine for stars with 
mass larger than 1.8 M0. 

The discretization of the PDE is based on a spectral decomposition: Cheby- 
shev polynomials radially and spherical harmonics horizontally. The distorted 
shape of the star is managed through a change of variables mapping the star to 
the spherical geometry. Internal precision is monitored through spectral conver- 



^The ESTER code is a public domain code that is freely available at 
:p: //code .google . com/p/ester-project/ _ 
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Table 1. Comparison between observationally derived parameters of the stars a Oph, a 
Lyr and a Leo and ESTER models. Data are respectively from iMonnier et"aD l|2010l l. 
iMonnier et aP \2'dVt \ and lChe et all l|201ll ). Note the good matching of the models. 



Star Ras Alhague (a Oph) Vega (a Lyr) Regulus (a Leo) 

Obs. Model Obs. Model Obs. Model 



Mass (M0) 


/I+0.23 
^•^-0.37 


2.22 


9 T c+0.10 


2.374 


4.15±0.06 


4.10 


Req (R©) 


2.858±0.015 


2.865 


2.726±0.006 


2.726 


4.21±0.07 


4.24 


Rpol (R-o) 


2.388±0.013 


2.385 


2.418±0.012 


2.418 


3.22±0.05 


3.23 


Tcq (K) 


7570±124 


7674 


8910±130 


8973 


11010±520 


11175 


Tpoi(K) 


9384±154 


9236 


10070±90 


10070 


14520±690 


14567 


L(L0) 


31.3±1 


31.1 


47.2±2 


48.0 


341±27 


351 


Veq (km/s) 


240±12 


242 


197±23 


205 


336±24 


335 


Peq (days) 




0.598 




0.672 




0.641 


Ppol (days) 




0.616 




0.697 




0.658 






0.70 




0.7546 




0.70 


-^corc / -^cnv . 




0.37 




0.271 




0.5 


z 




0.02 




0.0093 




0.02 



gence, the virial test and the energy test (see iRieutord fc Espinosa Larall2012l . for 
details). Iterations follow the Newton algorithm. Fig. [T] illustrates the result of 
a two-dimensional model for a massive star rotating at 98% of the critical angu- 
lar velocity. We show here the differential rotation as a function of radius and 
latitude. 



6.2 Gravity darkening law 

One of the first results of steady models has been the prediction of the gravity 
darkening law. Until now, the standard recipe was the von Zeipel law stating that 
the effective temperature is proportional to the j power of the effective gravity. 
Fitting brightness distributions, interferometric data of rapidly rotating stars have 
shown that other laws like T^s oc g^g with adjusted /3, were more appropriate. 
This has also been the conclusion of ESTER models with predictions on the /? 
value, although th e models show that a power law i s not exactly representing the 
gravity darkening ( Espinosa Lara fc Rieutord 2011 , 2012f) . In Fig. [51 we show the 



dependence of the flux with colatitude for a 30 M© star rotating close to critical 
angular velocity. The von Zeipel law underestimate the flux in the equatorial 
regions by more than a factor 3. 



6.3 Some models of nearby fast rotating stars 

Fundamentals parameters of some nearby (famous) rotating stars, derived by inter- 
ferometry, have been compared successfully to ESTER models as shown by Tab.fH 
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Table 2. Comparison between observationally derived parameters of the stars and ten- 
tative two-dimensional models. Data from 5 Vel are from iMerand et aD (|201lh . those of 
Achernar are from iDomiciano de Souza et~al] (|2012l ). The models compare nicely with 
observationally constrained data for the two components of S Vel A (an eclipsing binary) 
but have difhculties with Achernar. 



Star Delta Velorum Aa Delta Velorum Ab Achernar {a Eri) 

Obs. Model Obs. Model Obs. Model 



Mass (Mq) 


2.43 ±0.02 


2.43 


2.27±0.02 


2.27 




8.20 


Rcq (R-o) 


2.97±0.02 


2.95 


2.52±0.03 


2.52 


11.6±0.3 


11.5 


Rpol (R-o) 


2.79±0.04 


2.77 


2.37±0.02 


2.36 


8.0±0.4 


7.9 


Teq (K) 


9450 


9440 


9560 


9477 


9955_233g 


11250 


Tpol(K) 


10100 


10044 


10120 


10115 


18013^^:^^ 


16800 


L (Lq) 


67±3 


65.2 


51±2 


48.5 


4500±300 


3700 


Vcq (km/s) 


143 


143 


150 


153 


298±9 


339 


Peq (days) 




1.045 




0.832 




1.72 


Ppoi (days) 




1.084 




0.924 




1.68 


-^env. 




0.70 




0.70 




0.74 


-^corc/-^cnv. 




0.10 




0.30 




0.05 


z 




0.011 




0.011 




0.04 



A further comparison is shown in Tab. [2l S Vel A is a.n eclipsing binary that ha s 
been studied in detail by interferometry ([Merand et al.lbOlltlPribulla et al.ll201ll ). 
Here too, two dimensional models nicely fit the observationnally derived param- 
eters of the two startH. In addition they suggest that they are less metallic than 
the Sun (Z~ 0.011). We also note the stronger hydrogen depletion in the core of 
the most massive one as expected from its fastest evolution. 

On the other hand the case of Achernar, the closest Be star to the solar system, 
is more difficult since none of the explored models really fit the parameters derived 
from interferometry. Fitting the polar and equatorial radii lead to rather extreme 
composition or mass suggesting that the star may just have left the main sequence 
(to which we are c onstrained at the mo ment). In view of previous results based on 



spherical models (|Vinicius et al.l 120061 ) . this not totally surprising. However, this 



star remains a challenging case for two-dimensional models. 



7 Outlooks 



Presently, ESTER two dimensional models are, strictly speaking, models of in- 
ternal structure of rotating stars. No time evolution is included. They describe 
main sequence early- type stars, that is for masses larger than 1.8 Mq. Evolution 



■^The tidal distortion is weak, less than 10 
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along the main sequence can be mimicked by varying the hydrogen mass fraction 
in the convective core. Clearly, the next important steps are the extension to low 
mass stars, which means the computation of outer convective envelope, and the 
inclusion of time evolution. Other t wo-dimensional models are current ly being 
developed by Deupree and coworkers (|Deupreell201ll : iDeupree et al.ll2012l) . but in 
these models the differential rotation should be prescribed. 

Such internal structure models are useful to interpret the seismic frequencies 
of rotating stars when perturbative methods f ail. Presently, two code s may deal 
with two dim ensional models: TOP by Reese (jReese et al.ll2006L 120131 ) or ACOR 



by Ouazzani (jOuazzani et al.ll2Q12l ). Beside asterosismology, 2D- models are im- 



portant for interferometry. Indeed, the interpretation of interferometric visibilities 
requires the adjustement of models that include the c entrifugal distortion of the 



stars as well as the associated gravity darkening (e.g. iDomiciano de Souza et al 
2002[ ). Here, a future improvement will be the calculation of two-dimensional 
models of atmospheres. 

Beyond the obvious improvements that have been mentioned above, we see 
that progresses in the understanding of rotating stars will have to go through a 
better modeling of the turbulent transport in all places where turbulence develops. 
This is certainly the most challenging issue of this modeling since we do not have 
a general theory of turbulence at hands. 
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